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We present a statistical analysis of the two-point vorticity probability density of the vorticity field 
generated in the inverse cascade of stationary two-dimensional turbulence. 
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I. INTRODUCTION 

Fully developed turbulent flows are systems far from 
equilibrium giving rise to a transport of energy or enstro- 
phy across scales. Up to now, no generally accepted the- 
oretical description of this transport process has emerged 
although the phenomenological theory based on the 
works of Kolmogorov, Onsager, Heisenberg and their suc- 
cessors is quite successful in describing the gross features 
of turbulent fields (we refer to the monographs PQ, [2])- 

Two-dimensional turbulent flows play a prominent role 
due to the existence of a direct as well as an inverse 
cascade, as has been emphasized in the seminal work of 
Kraichnan [3] . The coexistence of both cascades has been 
convincingly demonstrated by Boffetta and Musacchio 
[I]. Recently, a detailed analysis of the contour lines of 
zero vorticity has revealed interesting scaling behaviour 
pointing to the existence of nontrivial multi-point statis- 
tics of vorticity in the inverse cascade [5] . The inverse cas- 
cade is a central issue in classical nonequilibrium physics 
0- 

The present contribution is concerned with the two- 
point statistics of vorticity in the inverse cascade. Start- 
ing from the hierarchy of evolution equations for multi- 
point probability distributions of vorticity, we shall per- 
form a detailed examination of the partial differential 
equation defining the probability distribution of the vor- 
ticity increment relying on input from direct numerical 
simulation (DNS). 



II. RESULTS FROM DNS 

The spatio-temporal dynamics of the turbulent field is 
governed by the two-dimensional vorticity equation 

d 

— w(x, t) + u(x, t) ■ V x w(x, t) = L(-AWx, t) + F(x, t) 
at 

(1) 

The velocity field u(x, t) is determined from vorticity 
w(x, t) via Biot-Savart's law, 




FIG. 1: Spectral energy flux obtained in DNS of 2d turbu- 
lence. The compensated spectrum exhibits a clear indication 
for the existence of the inverse energy cascade. The red line 
in the inset marks the band of Kolmogorov constants found 
in earlier publications [THS]. 



Viscosity is taken into account by the operator L(— A), 
where we allow for the presence of hyperviscosity as well 
as friction acting on the large scales. In order to achieve 
a constant energy flux in the inertial range we have used 
a hyperviscosity of the form vh{— A) 8 and a large scale 
friction — iv(— A) -1 to extract energy transported to the 
large scales by the inverse cascade. Energy is injected 
into the system by the forcing term F(x,i). Here we 
employ a forcing that acts only on a narrow wavenumber 
band around fco in Fourier space. Within this band the 
amplitudes |wk| of the Fourier modes 



u k (t) 



*k(t) 



(3) 



f x - x' 

u(x,i) = / dx'e^ x — -^w(x',t) 

J Z7T|X — x'l^ 
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are kept constant by resetting them every time step to the 
constant wo while the phases can evolve freely according 
to the vorticity equation [8] . In Fourier-space, the evolu- 
tion equations can be written for amplitudes and phases 
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as 



(Q(r) 2 ) 1/2 H(Q,r) 



{Cl{rf) ll% lel 



dt 



i . d . 



W k = L(-fc 2 )K| 
^(fc)-K|]6(fc,fc ) 
M k = 



(4) 



Here, N^, M k are the real and imaginary parts of the 
Fourier transform of the convective term u • Vw. F(k) 
is a real function centered around fco and T — > is a 
time relaxation constant. The function 0(fc, fc ) is dif- 
ferent from zero only in the forcing interval. From this 
representation it becomes obvious that the forcing can 
be considered to be deterministic. 

In most simulations (see e.g. [7HH]) the forcing is lo- 
cated at a wavenumber fco that is close to the highest 
wavenumber resolved . In combination with hyperviscos- 
ity this allows for a wide inertial range and suppression 
of the direct enstrophy cascade. However, this leads to 
a poor numerical resolution of the forcing scale. To cir- 
cumvent this flaw we inject the energy at an intermediate 
length scale and used a hyperviscosity that, for numer- 
ical reasons saturates at a maximal value. Comparison 
with other numerical simulations showed that this pro- 
cedure does not affect the results while yielding a better 
resolution at the forcing scale. 

The vorticity equation is numerically solved by means 
of a standard pseudospectral method with 2/3-dealiasing 
in a periodic box with side length 2n and a resolution 
of 2048 2 grid points. Time stepping is performed by a 
memory saving Runge-Kutta scheme of third order |10j . 
The hyperviscous term is treated by an integrating factor. 

Data-analysis has been performed over more than 200 
snapshots of the vorticity field taken from a simulation in 
a statistically stationary state. The snapshots are equally 
separated by about 1.4 large-eddy turnover times. 

Figure ([!]) exhibits the energy spectrum in Fourier- 
space in the inverse cascade, as obtained from our nu- 
merical solution. The spectral energy flux is constant and 
demonstrates the existence of the inverse cascade. The 
obtained Kolmogorov constant K a , defined by E(k) = 
K e 2 / 3 k~ 5 / 3 is consistent with the results reported in 
the literature, see e.g. [THS]- The probability distribu- 
tions of the vorticity increments, defined according to 
O x (r, t) — u)(x + r, t) — oj(x, t), are exhibited in fig. ([2| 
as a function of scale r — |r|. Apparently, the vorticity 
increment is nearly Gaussian distributed. Furthermore, 
the second moment of the vorticity increment exhibits 
oscillations as a function of the distance r, which are 
clearly a result of the small but visible oscillations in the 
shape of the pdfs. It may well be that the wings of the 
pdfs are super-Gaussian. In the following we shall be 
concerned with the behaviour represented by the center 
of the pdf, where a Gaussian representation of the statis- 
tics is possible. The behaviour of the tails of the pdf 
may be accessible by instanton calculations, as has been 
discussed in, e.g., [TT] . 



10 




1 f\"i 

l.uz 




■ 1 1 1 1 1 1 1 - 




1 

10" 1 


/ \ 


1.01 


-i 


10" 2 




1.00 
0.99 




10" 3 




0.98 




10" 4 


1 ^\ 


0.97 




10" 5 




0.96 
0.95 




10 -6 




0.94 






-6 -4 -2 2 4 6 



FIG. 2: Left: Probability distribution of the scaled vor- 
ticity increment f2(r, t)/v/ t) 2 ) as a function of r 
(r = {1, 2, . . . , 5} A). The Taylor scale A is defined by 
[2(u 2 )/({d x u) 2 )] 1/2 . Right: Second moment (tl(r,t) 2 ) of the 
vorticity increment as a function of scale r. The moment 
exhibits pronounced oscillations in scale. The constant ctq 3 
denotes the value of the second moment in the limit r — > oo, 
i.e. a£ = (2(c 2 )) 1 /2 



An important statistical relationship is the analog of 
the von Karman-Howarth equation formulated for vor- 
ticity. It explicitly reads 

V Xl • (u(x 1 ,t)u(x 1 ,t)u(x 2 ,t)) + 
V X2 • (u(x 2 ,t)uj(x 1 ,t)uj(x 2 ,t)) = 
(u(x l ,t)L(-A^ 2 )u}(x 2 ,t)) 
(w(x a ,t)L(-A Xl Mxi,t)) 



(u(x u t)F(x 2 ,t)) + (u(x 2 ,t)F(x u t)) 



(5) 



and can be derived in a straightforward manner from 
the evolution equation for vorticity. The transport terms 
of the left hand side of the equation are responsible for 
the transport process in the cascade. They essentially 
depend on the three-point vorticity correlations 



(w(x',t)w(xi,t)w(x 2 ,i)) 



(6) 



since the velocity field is determined from Biot-Savart's 
law. 



III. STATISTICAL TREATMENT 



We shall base our statistical description on the N- 
point vorticity pdf f N ({uii, Xj}, t), which we define as an 
ensemble average (f N ({uJi, x^}, t)) over the fine-grained 
probability distribution / 



xj, t) = nf =1 S(^ - tufa, t)) 



(7) 
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Here, w(xj,f) denotes the solution of the vorticity equa- 
tion at spatial point Xj and time t. 

It is straightforward to derive the following evolution 
equation for the fine-grained pdf (see e.g. [H]). A short 
discussion of the procedure can be found in the appendix. 
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N 

Z— 1 



/({WjjXi},*) 
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[L(-A x >(xj) + F(xi,t)] f({ui, Xj}, t) (8) 



In order to obtain the evolution equation for the coarse- 
grained pdf f N ({u>i,Xi},t), one has to perform an en- 
semble average of this equation. This averaging proce- 
dure immediately confronts us with the closure problem 
of turbulence, since the following unclosed expectations 
arise 

(u(x,t)f({u)i,Xi},t)) 

= {u,x\{Ui,Xi},t)f({uJi,Xi},t) 

(L(-A x )u(x,i)/({^,xJ,t)) 

= (£(- A ) U : x|{Wj, X 4 }, t)/({cjj, Xj}, t) 

(F(x,t)f({cJi,Xi},t)) 

= (F,x\{u i ,x i },t)f({cj i ,x i },t) (9) 

It is convenient to express these unclosed expectations, as 
shown above, in terms of conditional averages of velocity, 
(u, x\{u)i,Xi},t), dissipation (L(— A)u, x|{wj, Xj}, t) and 
forcing (F, x|{wj, Xj}, t) . The introduction of these quan- 
tities is justified by the fact that conditional averages 
can be estimated from DNS (see e.g. [13HT6] . and for a 
recent application to turbulent Rayleigh-Benard convec- 
tion [T7j). For the case of a white noise forcing, the term 
related to the conditional forcing may contain second- 
order derivatives with respect to u>i and then takes the 
form of a diffusion operator, Y^ij ^jrQi^i ~ x j) g£r/- An 
alternative way is to write down the whole chain of evo- 
lution equations for higher N-point functions, a proce- 
dure which is well documented in the literature [Tl [TBH2"0"] . 
However, up to now conclusive truncations of the hierar- 
chy are lacking. Additionally, we want to mention that a 
similar approach for the study of the velocity increment 
statistics of the inverse cascade has been undertaken by 
Boffetta et al. [2T] based on theoretical approaches by 
V. Yakhot [22], [23] . Here, numerical input for the condi- 
tional pressure term is used to close the equation for the 
probability distribution of the velocity increment. 

Using the conditional expectations ^ we arrive at the 
following partial differential equation determining the pdf 



d_ 
at 



N 



£ V *» • (u,Xi\{u)i,Xi},t) 



i=l 
N 



f({Ui,Xi},t) 



i=l 

+ (i ;l ,x 4 |{a; 4 ,x. i }),i]/({w i ,x l },t) 



(10) 



These partial differential equations can be solved by the 
method of characteristics. The corresponding character- 
istic equations read 

Xj = (UjXil-fu^Xi},*) 

ibi = (i(-A)w,Xj|{a;j,x^},t) + (F,Xi\{ui,in},t) (11) 

Along these characteristic curves the probability distri- 
bution changes according to 



dt 



f({uJi,Xi},t) 



Y V Xj • (u,Xi\{u h Xi},t) 

i 

d 



/({Wj,Xj},t) 



£ 

i 

£ 



_d_ 



(L(-A)u),Xi\{u)i,Xi},t) 



f({Ui,Xi},t) 



(FjXjljcjj.Xj},*) 



f({(Ji, Xj}, t) 







(12) 



For the case of a stochastic forcing, the characteristic 
equations are stochastic differential equations. In the 
following we shall also include a discussion of this case, 
where the characteristic equations take the form 



Xj 



(u,Xj|{u;j,Xj},t) 

(L(-A)u,x i \{u i ,x i },t} + F i (t) (13) 



Here, Fi(t) are stochastic white noise forces with corre- 
lation function 



(F l (t)F J {t r )) = 2Q(xj - Xj)6{t - t') 



(14) 



In fact, the expression for the conditional force will allow 
such an interpretation. In the following we shall discuss 
an interpretation of these stochastic differential equations 
as a kind of stochastic point vortex model. 



IV. LINEAR GAUSSIAN APPROXIMATION 

Since the two-point statistics of vorticity is close to 
Gaussian, it is tempting to approximate the conditional 
averages Q based on the assumption of Gaussian statis- 
tics. The conditional velocity field can be expressed via 
Biot-Savart's law in terms of the conditional vorticity ex- 
pectation (oj, x'|{u;j,Xj},t) according to 



(i^xl-ju^Xj},*) 
dx'e, x 



-( W ,x'|{ Wl , Xl },t) (15) 



2ir x - x : 



In the case of Gaussian statistics all quantities depend 
on the vorticity-vorticity correlation function which we 
denote by C(x — x'). However, for Gaussian statistics 
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FIG. 3: Conditional vorticity and velocity fields, 
(id, x\uii, asi; Wa, 052) and (it, £c|cji, cci; CJ2, 232} respectively, for 
the case wi = — u>2 ~ (oJ 2 ) 1 ^ 2 and r = |a52 — xi| = X obtained 
from DNS. The vorticity is color-coded, while the velocity is 
represented by arrows and its streamlines. 



the conditional expectation is given by (c.f. [12] and the 
appendix of this paper) 



x\{u}j,Xi},t) =y^(7(x-x fe )(7(x fc -x;) l u>i (16) 



so that we obtain the following linear approximations for 
conditional velocity and conditional vorticity dissipation 

(u, x|{wi, xj, t) = £ U(x - x fc )C(x fc - x,)" 1 ^ (17) 

kl 



(L(-A)w, x|{wj, xj, i) = -y 7(x - x fe )C(x fe - xj) x u;; 

(18) 



;,7 



We have introduced the velocity kernel 
U(x-x,) = ^e 2 x/dx'^C(x'- Xl) 



^(Ix-Xil) 



(19) 



The field U(x — Xj) is the velocity field of a vorticity 
distribution C(x — x.j) centered at Xj. A point vortex 
with circulation Ti would be described by C(x — Xj) = 
TiS(x — Xj). As a consequence we can visualize the ve- 
locity field U(x — Xj) as a field of a screened vortex, i.e. 
a quasi-vortex in the sense of a Landau quasi-particle 
taking into account the surrounding turbulence by an 
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FIG. 4: Conditional vorticity and velocity fields as in Fig. [3] 
but here evaluated following the linear Gaussian approxima- 
tion, Eqs. ( 16 1 and ( 17 1 respectively. 



effective velocity profile. The complete conditional veloc- 
ity field is a linear superposition of velocity fields of such 
screened (or dressed) vortices. Alternatively, one could 
also say that Biot-Savart's law is changed by the effec- 
tive consideration of the turbulent background field. We 
mention that a change of Biot-Savart's law has been used 
by Chevillard et al. [37] in the modeling of the statistics 
of velocity increments in three dimensional turbulence. 

Let us explicitly determine the conditional velocity 
field conditioned on vorticity uj\, L02 at positions Xi, X2: 



(u,x|{wi,: 
= U(x- 



i}) 



C(0)cji - C(xi - x 2 )w 2 



+u(x - x 2 ; 



C(0) 2 
C(0)u 2 - 



-C( Xl 
C(x 2 - 



-x 2 )2 

Xi)wi 



(20) 



C(0) 2 -C( Xl -x 2 ) 2 
Furthermore, we note that U(0) = 0, due to the fact 



that C(0) is finite, c.f. eq. ([19 1. 

Furthermore, a dissipation field 7(x 
according to 



Xj) is defined 



7 (x- Xi ) = -L(-A)C(x- Xi ) (21) 
The characteristic equations then take the form 
x 4 = ^ U ( x * - x fe )C(x fe - X() _1 a; ; 

kl 

lo 1 = -Y j - x fe )C(x fe - x,) -1 ^ + F( Xi ,t) (22) 



We have determined the conditional vorticity field from 
DNS (for the case of two points) and have compared it 
with the one obtained from the linear Gaussian approx- 
imation. Visually, there is no difference, as can be seen 
from Figs. Q and d6j), respectively. However, 
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FIG. 5: Same quantities as in Fig.[3] now for r ~ 4A. To 
be able to show the conditional averages for large r in terms 
of the Taylor scale which is still deep in the inertial range 
we used a 2048 2 simulation, where the forcing scale is less 
resolved. 



as we shall see, the linear Gaussian approximation is not 
able to describe the turbulent cascade, at least if one 
restricts oneself to the treatment of two-point statistics. 
The reason is that the characteristic equation for the case 
N = 2 takes the form 

Xl = U(Xi - X 2 )G(W1,W2, |xi - x 2 |) 

x 2 = U(x 2 - xi)G(w 2 ,wi, |xi - x 2 |) (23) 

with 

nt i w G(0)cj 2 - G( Xl - XgM 

G(c 1 ,o; 2 ,|x 1 -x 2 |)= c(0)2 _ c(xi _ X2)2 (24) 



It describes two-vortex motions of vortices in a 
screened velocity field. This motion does not allow for 
a variation of the distance, i.e. there is no transport pro- 
cess in scale. This can be seen by invoking the analog 
of the van Karman-Howard equation for vorticity. The 
transport terms read under consideration of the condi- 
tional velocity field 

J du x duj2 wio; 2 [U(xi - x 2 )G(w 2 , w t , |x x - x 2 |) • V Xl + 

U(x 2 - xi)G(w 2 , oji, |xi - x 2 |) • V X2 ]/(cji,Xi;w 2 ,x 2 ) 

(25) 

However, since for homogeneous isotropic statistics the 
probability distribution depends on |xi — x 2 |, the trans- 
port term vanishes since 
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FIG. 6: Linear Gaussian approximation as in Fig.|4]for r « 4A. 
The estimation was based upon the same simulation as in 
Fig.0 

and the gradient of f(u>i, Xi; w 2 , x 2 ) points in the direc- 
tion of the vector Xi — x 2 such that 



U(xi - x 2 ) • V Xl /(wi,Xi;a; 2 ,x 2 ) = 



(27) 



The characteristic equations ( 22 ) can be interpreted in 



a Lagrangian sense. Neglecting the temporal variation 
of the vorticities tOi the set of differential equations for 
the positions Xj can be viewed as a modified point vortex 
system, where the velocity field of an individual vortex 
is not the one of a point vortex, but is the velocity field 
of a dressed vortex. This perspective can be seen in close 
analogy to the notion of a Landau quasi-particle. This 
notion is widely used in theories of many particle sys- 
tems, where a part of the many-particle interaction can 
be absorbed into effective properties of single particles. 
In the present case, the reduction from a continuum de- 
scription of vorticity to the description of the statistics of 
the field at finite spatial points quite naturally leads to an 
effective point vortex dynamics. However, since we are 
still in the realm of the linear Gaussian approximation, 
we are not able to obtain the inverse cascade on the ba- 
sis of a two-point closure due to the vanishing transport 
terms discussed above. One might think of obtaining the 
inverse cascade starting from the three vortex problem 
in the linear Gaussian approximation. The three-vortex 
problem certainly brings in new scales, however, we con- 
jecture that this is not the correct way to proceed. 

Due to the temporal variation of the vorticities u)i the 
characteristic equations form an extension of the problem 
of point vortex motion. We refer the reader to the review 
[2"1] . On the other side, if the vorticities are assumed to 
undergo white noise processes 



(uii^ujjit')) = 2Q ij 8{t-t') 



(28) 



U(xi - x 2 ) = e 2 x 



xi - x 2 
h(\xi - x 2 |) 



(26) 



one obtains an extension of the Kraichnan model |25j . 
which has been investigated with respect to passive scalar 
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advection (for a review we refer the reader to [25] ) . As is 
transparent from ( |22[ ), in the adiabatic limit the vortici- 
ties are directly related to the white noise forces 



10 



(29) 



where Tu is the inverse matrix to J2u 7( x * ~ x fc)C( x fc — 
x;) _1 . The result is the Fokker-Planck equation for the 
positions x; of the point vortices 



0_ 

at 



p({xi},t) = Y,V* i D ij [{x l }}Vx jP ({x i },t) (30) 



In general, the diffusion matrix Dij[{x;}] depends on all 
positions {x;} of the vortices. In the Kraichnan model, 
the diffusion matrix Dij is only a function of the differ- 
ence vector x.j — Xj. The explicit expression of fl y [{x[}] 
is 

Ai [{*}] = 

U(x, ; - x i )C i - 1 r mn U(x j - * k )C^T m Q m (31) 
We sum over dummy indices. 



V. VORTICITY INCREMENTS 

For the following we shall restrict our considerations to 
the case of two spatial points, since then the estimation 
of the conditional averages from DNS becomes feasible. 
Due to homogeneity and isotropy, the two-point pdf only 
depends on the distance r, r = xi — x 2 and the evolution 
equation for the two-point functions takes the form 



- — n{r,u)i,u) 2 ) + ^ — fJ.[r, w 2 ,wi) 

OUJi OUJ2 



(32) 



Thereby, we have defined the conditional velocity incre- 
ment 

v(r,wi,w 2 ) = (u,Xi | {u)i,Xi}) - (u,x 2 |{wi,xj) (33) 

and have taken into account that the quantity \i = 
L + F only depends on the distance r between the two 
points. Furthermore, due to radial symmetry of the pdf 
/(wi, u>2, r), only the radial component of the conditional 
velocity v r {u\, w 2 , r)e r contributes (e r denotes the unit 
vector in radial direction). It is convenient to introduce 
the pdf /i(wi,w 2 ,r) = rf{u\,U2,r). 

It is further possible to reduce the complexity of the 
problem by considering the statistics of the vorticity in- 
crement SI = cji — cj 2 . The corresponding equation defin- 
ing the statistics of Q, i?(Q,r) can be obtained via the 
definition 



H(Q,r) = /dwidw 2 £(ft-(wi-W2))Mwi,W2,r) (34) 
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FIG. 7: The function g(r) and the fit according to eq. (38 1. 
Scale is measured in terms of the Taylor scale A. The integral 
scale is denoted by L. As above we used the constants <rj° = 
(2(u 2 )) 1 ^ 2 and <jq = (2(oj 2 )) 1//2 respectively to normalize our 
results. 



A straightforward calculation yields 



d 
Or 



V r (r, Q)H(Q,r) 



d_ 

an 



/i(r,n)H{n,r) (35) 



The conditional velocity field V r (r, Q) has to be an even 
function of SI, furthermore its average has to vanish: 

V r (r,Q) = -V r (r,-fi) 

Jdfl H(fl,r)Vr{r,Q) = (36) 

We have determined the conditional velocity field 
V r (VL, r ) as well as the conditional field /i(Sl, r) from DNS. 
For not too small values of r it is possible to perform the 
following quadratic fit in SI 2 , see fig. (|7|): 



V r (Q,r) = g(r)[fl 2 - (SIM) 2 )] 
ju(0, r) = fj,(r)Cl 



(37) 



This fit is consistent with the conditions ( 36 1 and can be 
viewed as a Taylor expansion in Q. The quantity /Li(f2, r) 
can be represented by a linear fit in SI. 

It turns out that the function g(r) is an oscillating func- 
tion of r (see fig. (j7|)) which can be well approximated in 
the inertial range by 



9{r) 



' sin (At — $) 



(38) 



The constant k = 2.1-k/X is related to the width of the 
forcing function of the turbulence. The exponent 8 is of 
the order of 5 ~ 0.8. Amplitude and phase turn out to 
be c = 1.8 x IOV^/ct^ 2 , $ = 1.3vr. 



The function /x(r, SI) consists of three contributions, 
namely forcing, dissipation, and friction. All contribute 
with a linear dependency in SI to /i(r, SI). The various 
contributions are exhibited in fig. 
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FIG. 8: The three contributions to the function fj,(r, Q), fitted 
by n(r,£l) = n(r)Q, due to forcing, viscosity, and friction. 



FIG. 9: The function fi(r) for the contributions due to 
forcing, viscosity, and friction compared with the relation 
fi(r) — {Q(r) 2 )g' (r). Scale is measured in terms of the Taylor 
scale. 



We also have performed a fit of the form 



C(r) = ar a [c + sin(W + $)] 



(39) 



with the parameters a — 0.0016, to — 104, c = 0.66, 
$ = 0.32 and a = —4/3 to the correlation function 
explicitly incorporating the expectation of scaling be- 
haviour according to K41. 

It is interesting to notice that we can find a solution 
to the pdf equation ( 35 ) with the ansatz ( 37 ) for values 
of the scale, for which the function g(r) vanishes. In 
that case the characteristic equation degenerates to the 
differential equation 



/(r)[n 2 -<^ 2 (r,*)>]#aV) 



d_ 

"on 



li(r)nH(n,r) (40) 



with the explicit Gaussian solution 
H(n,r) 



e 2<Q(r.t)2> 



(41) 



The following relationship between conditional velocity 
field, conditional dissipation and forcing, and the corre- 
lation function (fi(r, t) 2 ) then has to hold: 



(42) 



A check of this equality is presented in fig. We want 
to point out that, strictly, the result is only valid for 
values of r, for which g(r) = 0, i.e. at the minima and 
maxima of the moment. In between these points, the pdf 
can be determined perturbatively. 




0.0 10 2.0 3.0 40 50 6.0 



FIG. 10: Characteristic curves of the characteristic equation 
(221. The vertical lines mark the scale with g(r) — 0. The 
curve is the correlation (f2(r, t) 2 ). 



VI. LAGRANGIAN MODEL 

It is plausible to look for a relationship of the char- 
acteristic equation (111 with the Lagrangian description 



of the relative motion of two-point particles with relative 
vorticity f2. The conditional velocity field thereby is a 
type of averaged relative velocity field, whose single re- 
alization has a deterministic as well as a stochastic part. 
In order to formulate a model for these single stochas- 
tic realizations, we separate the conditional velocity field 
into two parts, a deterministic and a fluctuating part: 



V r (r,Q) = G{r)[n 2 
- D(r) 



H{r, fi) dr 



d 

H(r,n) 



(43) 



If we take the Gaussian solution (41) valid for the scale 
separation defined by g(r) = 0, we are led to the relation 

1 d 



g(r) = G(r)-D(r) 



G(r) 



2(n(r,t) 2 ) 2 dr 
D(r) 1 



(Q(r,ty 
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(44) 



At the present stage, we are not able to fix the function 
D(r). Therefore, we take D(r) as a small positive quan- 
tity in the following. A plausible choice would be to take 
D(r) according to Kraichnan's model of passive scalar 
advection 2Bj. 

A similar ansatz is made for the quantity /Lt(r, f2), i.e. 
we assume a linear damping and white noise forcing and 
model the Lagrangian statistics by the Langevin equa- 
tions for relative position r and relative vorticity f2: 

r = G(r)[n 2 -(Q(r,t) 2 )]+r, 

fi = -7(r)fi + F (45) 



We intend to derive the decomposition (44 1 and the 



statistics of F using well-established methods for the de- 
termination of drift- and diffusion terms for diffusion pro- 
cesses [28 directly from Lagrangian data. 

Thereby, the two stochastic forces rj, F are white 
noise forces and are characterized by their correla- 
tions (F(r,t)F(r,t')) = Q(r)S(t - t'), ( V (r,t)r,(r,t')) = 
C(r)8(t — tf). A comparison yields the relationship 



fx(r) = -7(r) + 



Q(r) 



(46) 



(Q(r,*) 2 ) 

We have used the Ito-interpretation. The correspond- 
ing Fokker-Planck equation for the Lagrangian joint 
position- vorticity increment H(r,£l,t) reads 

d d 



dt dr 



d_ 

on 



7(r)n + Q(r 



d 2 

on 2 



H(r,n,t) 



H(r,Sl,t) 



(47) 



It is interesting to note that the Langevin equations are 
prototype equations for so-called ratchets. The ratchet 
process is a genuinely nonlinear transport process under- 
lying the motion of Brownian motors and rotors (we refer 
the reader to the reviews [3D]). According to [3"T] . 
the present system is a pulsating ratchet. 



It is interesting to compare the Lagrangian model ( 47 ) 
with the models presented by Castiglione and Pumir [34] 
and Pumir et al. |35j . who are dealing with multi-point 
statistics. Our present analysis refines these models by 
an explicit consideration of the functional form of the 
conditional relative velocity field between two Lagrangian 
points. 

On a more general level it would be interesting to un- 
derstand the statistics of Lagrangian observables like ve- 
locity increments and acceleration, which have been in- 
vestigated recently in two-dimensional flows [3H UDJ 
|4"T] in terms of the multi-point vorticity statistics. 



VII. SUMMARY 

A generalization of the Langevin Model to the two- 
point vorticity equation is straightforward and has the 
following structure 

xi = e<pU(r, u)i,u}2) + e r V(r,ux,0J2) + m 

x 2 = -e v U(r,u>2,LJi) — e r V(r,wi,u)2) +»?2 (48) 

The vorticity variables undergo Ornstein-Uhlenbeck pro- 
cesses. U (r, u>\i 0J2) is an odd function with respect to uj.i , 
whereas V(r, 1^1,0:2) is an even function. In the present 
paper we have identified and used the lowest order ap- 
proximations 

e v C/(r,wi, lo 2 ) = U(xi - x 2 ) C^/wj 

i 

V(r,ux,U2) = G(r)22aij(r)[ujiU)j-(wiCjj}](A9) 

ij 

Our analysis has led us to an effective point vortex de- 
scription, which is in the spirit of the Landau effective 
particle approach for many-body systems. The effective 
azimuthal dynamics is a two vortex motion, however with 
a screened velocity field. The effective relative motion, 
which is absent for bare point vortices as well as for the 
vortex motion in the linear Gaussian approximation is 
the fingerprint of the inverse cascade. It determines the 
energy flux. The relative motion can be attracting or 
repelling, depending on the distance and the instanta- 
neous vorticities. A directed flux results due to a ratchet 
effect, a genuinely nonlinear transport process. It re- 
mains a challenge to derive this relative motion between 
two effective point vortices directly from the stochastic 
Navier-Stokes dynamics. 



Appendix A: Kinetic Equations 

For the sake of completeness, we present the derivation 
of the hierarchy of evolution equations for the vorticity 
probability distributions and explicitly discuss the intro- 
duction of the conditional averages. To this end we define 
the fine-grained N-point vorticity distributions 

/(wi,xi;i) = 8{ui - cj(xi,*)) 
/({wi.Xi};*) = Ili^Wi-wfo.t)) (Al) 

The evolution equations for these quantities are obtained 
by differentiation with respect to time: 

d d d 

^/(wi,xi;i) = --^cj(x l ,t)—S(ui 1 - uj(x 1 ,t)) (A2) 

Taking into account the vorticity equation we obtain 
d 



dt 



+ u(xi,t) • V X1 



d 
<9wi 



[£(-A)w(xi,t) + F(x 1> t)]/(wi J xi;t) (A3) 
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Thereby, we have used the fact that the velocity field is 
incompressible 



d 
duji 



u(xi, t) ■ V Xl w(xi, t)5(u>i - w(xi, i)) 



-u(xi,t) • V Xl /(wi,xi;i) 



(A4) 



We now proceed to derive the evolution equations for 
the coarse-grained probability distribution. Averaging 
yields 



0_ 



/(u>i,xi;f) +V Xl • (u(xi,t)/(wi,xi;i)) 



d 

d 
du}\ 



(F(x 1 ,t)/(o; 1) x 1 ;i)) 

^£(-A)w( Xl , t)/( Wl> x i; *)) (A5) 



The unclosed expectation value (u(xi,i)/(a;i,Xi;t)), 
can be expressed in terms of the conditional expectation 

(u(xi,i)/(cJi,Xi;t)) = U(xi|wixi;t)/(wi,xi;t) (A6) 

or, alternatively, using the representation 

u(xi, t) = J dx' U(xi - x')w(x', t) (A7) 



with 



U(x) = e z x 



(A8) 



2tt|x| 2 

in terms of the two-point probability distribution 
(u(xi,t)/(wi,xi;t)> = 

J dx' y*dw' w'U(xi -x'),/V,x';wi,Xi;i) = 

/" dx' /" dw' U(xi - x') (u, x '\ui , X! , t)/(«! , xi ; t) 

(A9) 

In a similar way, the dissipative term can be expressed 

as 

(i(-A a . 1 )w(x 1 ,t)/(«i,xi;t)) = 

ydx' Jdu>' 5( Xl -x')a/L(-A a /)/(w',x , ;« 1 ,x 1 ;t) 

(A10) 

Finally, we treat the forcing term. To this end we as- 
sume the force to be (^-correlated in time. The evolution 
equation takes the form 



This leads to 



(All) 



(A12) 



where we have defined 

Q(x 1 -x 3 -)*(t-f / ) = i(F( Xl ,t)F( Xj ,i)) (A13) 

This completes the derivation of the kinetic equation for 
the one-point probability distribution. 

The derivation of the equation for the multi-point dis- 
tribution is performed analogously. We obtain 



d_ 

at 



•V X4 ■V{x i \{oj h xi}) 
d 



f({wi,Xi};t) = 



i 

+(F,Xi\{ui h xi},t)}f({u)i,Xi};t) 
where the operator Li is defined according to 



Li = — ju>i 

OUli 



^Q(x 4 -Xj-) 



dujiUij 



(A14) 



(A15) 



We remind the reader that, in general, the conditional ve- 
locity field U(xi|{cL>z,xj}) does not need to be solenoidal. 



Appendix B: Conditional Expectations: Linear 
Gaussian Approximation 

It is straightforward to evaluate the expectation nec- 
essary to calculate the conditional velocity field based on 
the assumption of Gaussian statistics using the Fourier 
representation of the distribution [T2] 

/( w ',x',{^, Xi }H 

y da' n;da i e i ( aV+s; ' a " J, )e- w(Q: '' x '' {Qi ' Xi » (Bl) 

The quantity W(a' , x', {ai, x/}) is the cumulant generat- 
ing function of the multi-point vorticity statistics. The 
conditional expectation of vorticity is then obtained by 



du' w f(w ,x ;wi,xi; ..; ui N , x N ; t) 



Uidai 



]~W(a' ,x' ,{a hXl }) 



i oa i oloi 



/({««, x,}) (B2) 



This is an exact relationship involving the generating 
function W(a', x', {a/, x;}) of the (N + l)-point cumu- 
lants of vorticity. 

For a Gaussian distribution we have 



W(a',x',{a i ,x i }) = ^-C(0) 

+ ^^i c(x ._ x . )+ ^ a / a . c , (x /_ x . ) (B3) 

i,3 i 
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and we obtain for the expectation 

du/ w'/(u/,x'; {uji,xi}) = 

= J Hidai e jE " lQl x 

x 1 y C(x' - xi) ai e- S«J ^ c ^-^) 
i 

% 

= ]TC*(x' - xOC^C^})^/^^,^}) (B4) 



We have used the identity 



■E.,,^^-^) (B6) 



Here, CV, = C(xj — Xj) denotes the vorticity-vorticity 
correlation function at points Xj, Xj, which completely 
determines a multivariate Gaussian vorticity correlation. 
C^({x;}) denotes the inverse matrix, 



^C(x l -x fc )Cj l ({x ; })=4 



(B5) 



This expression for the conditional vorticity can be 
used to calculate the conditional velocity (17) as well as 



the conditional dissipation (21) in the Gaussian approx- 



imation. Since the conditional velocity field turns out to 
be linear in the vorticities, we shall denote the present 
approximation as linear Gaussian approximation (LGA) . 
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